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I.  INTRODUCTION 

Particle  size  diagnostics  are  Important  In  combustion  systems  that 
produce  and/or  consume  particles  and/or  droplets.  Production  of  soot  in 
the  burning  of  hydrocarbons  reduces  combustion  efficiency,  and  the  exhaust 
of  sooting  systems  is  a  source  of  air  pollution.  Measurement  of  soot 
production  in  flames  is  therefore  crucial  for  optimizing  the  efficiency  of 
combustion  systems  and  for  reducing  their  environmental  impact.  Data  that 
are  fundamental  to  optimizing  coal  and  liquid-fuel  combustion  systems  are 
the  ablation  rates  of  fuel  and  production  rates  of  soot.  Three  techniques 
for  obtaining  such  data  are  discussed  in  this  report. 

Optical  diagnostic  techniques  are  especially  suited  to  combustion 
systems,  because  chemically  reacting  flows  are  more  sensitive  to  the  pres¬ 
ence  of  invasive  probes  than  non-reacting  flows,  and  the  probes  themselves 
may  not  survive  in  the  hostile  combustion  environment.  The  mapping  of 
velocity  distributions  In  flames *»*  is  an  example  of  a  diagnostic  situa¬ 
tion  in  which  conventional  hot-wire  velocity  probes  cannot  be  readily  used 
but  In  which  use  of  the  laser-Doppler  velocimeter  (LDV)^  was  successful. 
Sampling  probes  for  subsequent  measurement  of  particle  size  In  sooting 
flames  are  suspected  to  modify  the  particle  size  and  even  its  shape. ^  In 
situ  light  scattering  techniques  avoid  such  problems  and  have  proven  to  be 
useful  as  combustion  diagnostics. 

Aerosol  particles  cover  a  size  range  of  at  least  four  orders  of  mag¬ 
nitude.  Particles  that  fall  within  the  first  two  orders  of  magnitude, 
between  10  nm  and  1  pm,  are  typically  smaller  than  the  wavelength  of  the 
interrogating  light.  Rayleigh  or  Mie  scattering  provides  a  signal  that 
contains  the  signature  of  the  scattering  particles.  These  smaller  par¬ 
ticles  undergo  significant  Brownian  motion,  which,  by  the  Doppler  effect, 
broadens  the  (initially  narrow)  spectrum  of  the  interrogating  light  as  it 
is  scattered  by  the  particle.  Measurement  of  this  spectrum  to  infer  par¬ 
ticle  size  is  known  as  diffusion-broadening  spectroscopy  (DBS).  DBS  is 
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reviewed  in  detail  in  Section  II.  Particles  in  the  second  two  orders  of 
magnitude,  between  1  and  100  pm  in  diameter,  are  often  amenable  to  photo¬ 
graphic  imaging  techniques.  Such  direct  measurements  of  particle  size  are 
preferred,  since,  unlike  DBS,  they  do  not  require  knowing  the  physical  and 
optical  parameters  of  the  particles.  Holographic  imaging  is  a  useful 
technique  in  this  size  range  and  is  discussed  in  Section  III. A.  A  final 
imaging  technique  developed  by  Chang  et  al . ^  uses  Raman  scattering  or 
laser  light  emitted  by  the  particle  itself,  and  is  described  in  Section 
III.B. 

Many  other  optical  techniques  are  useful  in  particle  sizing,  includ¬ 
ing  extinction  and  scattering  measurements  at  several  angles,  wavelengths, 
or  polarizations;  the  interruption  of  a  focused  beam;  and  LDV  measurement 
of  aerodynamic  particle  diameter.  This  report  emphasizes  DBS  measurement 
of  particle  size  because  it  is  a  point  measurement  that  can  be  accom¬ 
plished  in  situ  in  combustion  systems.  This  review  also  provides  a  com¬ 
plete  summary  of  the  experimental  and  theoretical  considerations  that 
apply  to  this  maturing  technology. 


II.  DIFFUSION -BROADENING  SPECTROSCOPY 


Soon  after  the  invention  of  the  laser,  Forrester^  discussed  the  pos¬ 
sibility  of  high-resolution  optical  spectroscopy  using  the  radar  technique 
of  beating  the  return  signal  with  itself  or  the  signal  of  a  local  oscilla¬ 
tor.  A  detector  that  responds  to  the  square  of  the  electric  field  (such 
as  a  photomultiplier)  yields  a  signal  whose  power  spectrum  contains  beats 
between  all  of  the  frequencies  present  in  the  optical  power  spectrum. 

This  technique  transfers  the  power  spectrum  of  the  scattered  light  from 
optical  frequencies  (10^  Hz)  to  the  0-10  MHz  range  and  is  the  basis  of 

both  LDV  and  DBS.  Many  other  applications  are  mentioned  in  review 
7—0 

articles. 

The  use  of  DBS  in  flames  is  reviewed  in  a  recent  article  by  Weill, 
Lhuissier,  and  Gouesbet*®  and  is  summarized  here.  DBS  was  first  used  to 
size  polymer  macromolecules  in  liquid  solutions,  for  which  Koppel 
proposed  the  method  of  cumulants  to  extract  information  about  the  particle 
size  distribution.  The  extension  of  DBS  to  aerosol  particles  was  accom¬ 
plished  by  Hinds  and  Reist,^  who  conducted  experiments  with  monodisperse 
aerosols  and  suggested  detectability  (the  signal-to-noise  ratio  at  0  Hz) 
measurements  for  polydisperse  systems.^  Edwards  et  al.^  discussed  the 
broadening  of  the  measured  spectrum  when  the  particles  are  in  a  flowing 
system. 

Penner,  Bernard,  and  Jerskey  2,16  appiie(j  DBS  to  soot  particle  sizing 
in  flames,  and  measured  particle  growth  on  a  premixed  C2H4-O2  flat 
flame. ^  Penner  and  Chang  present  further  experimental  results  and 
evidence  for  the  existence  of  self-preserving  particle  size  distributions 
in  flames,  concluding  that  detectability  measurements  are  not  sensitive 
enough  to  accurately  diagnose  log-normal  polydispersity,  except  for 
selected  size  ranges  using  polarized  beams.  Chang  and  Penner  also 
derived  the  effect  of  turbulence  on  the  measured  spectrum  and  on  the  abil¬ 
ity  to  measure  particle  size. 


Gouesbet  and  coworkers^®* ^  measured  mean  diameters  and  number  den¬ 
sities  in  CH4-O2  and  CgHg-C^  flat  flames  for  various  equivalence  ratios 
and  flow  rates,  comparing  DBS  with  particle  sizing  by  optical  extinction 
and  scattering  measurements*  From  detailed  maps  of  measured  number  den¬ 
sities  and  soot  volume  fractions,  they  concluded  that  coagulation  governs 
particle  growth  in  the  region  between  0.5  and  2.5  cm  above  the  burner. 

King,  Sorensen,  Lester,  and  Merklin  *  used  DBS  to  measure  mean 
sizes  and  number  densities,  and  proposed  the  method  of  cumulants  for  poly- 
disperse  systems  in  flames.  Taylor,  Scrivner,  Sorensen,  and  Merklin^ 
theoretically  compared  the  cumulant  method  with  inverse  LaPlace  transform 
techniques  for  diagnosing  log-normal  polydispersity  and  concluded  that 
both  methods  are  accurate  only  if  the  standard  deviation  of  the  particle 
size  distribution  is  less  than  1.4.  Recently,  those  authors  used  the 
cumulant  technique  and  measured  a  decrease  in  the  width  of  the  particle 
size  distribution  with  height  above  the  burner.  They  interpreted  the 
decrease  as  further  evidence  for  the  importance  of  coagulation  in  the 
early  regions  of  the  flame.  A  lean  flame  displayed  mostly  surface  growth 
after  1.5  cm.  Their  results  also  suggest  that  the  self-preserving  par¬ 
ticle  size  distribution  is  attained  in  these  flames. 

In  the  ten  years  since  its  first  use  in  flames,  DBS  has  been  demon¬ 
strated  to  yield  accurate  mean  particle  sizes  and  number  densities  for 
soot.  The  width  of  the  particle  size  distribution  has  also  been  measured 
through  detailed  analysis  of  the  photocurrent  correlation  function.  The 
limits  of  flow  velocity,  turbulence  level  and  scale,  and  size  distribution 
are  well  documented,  as  is  the  effect  of  extinction  and  scattering  on  the 
measured  parameters.  Uncertainties  in  the  flame  temperature,  viscosity, 
Knudsen  number,  and  particle  refractive  index  are  the  main  sources  of 
error,  but  the  ability  to  perform  measurements  in  situ  makes  DBS  the  pre¬ 
ferred  technique  for  measuring  soot  particle  size  in  flames. 

The  theoretical  considerations  on  which  the  above-described  experi¬ 
ments  are  based  are  detailed  in  Section  II. A. 


A.  THEORETICAL  BASIS  OF  DBS 

The  theory  of  DBS  is  based  on  the  broadening  of  a  laser's  optical 
spectrum  by  the  random  motion  of  many  diffusing  particles,  which  is  anal¬ 
ogous  to  the  diffusion  broadening  of  spectral  lines.  The  various  parti¬ 
cles  contribute  different  Doppler  shifts  to  the  light  they  scatter  in  the 
direction  of  the  detector,  and  the  duration  of  a  given  signal  is  inter¬ 
rupted  by  collisions  with  the  fluid  surrounding  the  particles.  The  time- 
averaged  detector  power  spectrum  yields  the  characteristic  broadening  of 
the  laser  light,  from  which  the  particle  size  can  be  inferred  by  means  of 
diffusion  theory. 

The  model  of  a  typical  laser-scattering  experiment  is  depicted  in 
Fig.  1,  which  also  shows  the  relation  between  the  wave  vectors  discussed 
below.  Only  homodyne  spectroscopy  (scattering  from  a  single  input  beam) 
is  considered  here.  Note  that  heterodyne  spectroscopy  (mixing  the 
horaodyne  signal  with  that  of  another  beam)  can  yield  information  about  the 
particle  velocity  (via  LDV)  as  well  as  its  size.^»1^ 

Consider  a  laser  beam  of  frequency  (i)q  propagating  in  the  direction  of 
the  wave  vector  Light  is  scattered  from  this  beam  (by  particles  pass¬ 

ing  through  it)  into  the  direction  of  a  photomultiplier,  denoted  by  £g. 

The  intersection  of  the  laser  beam  and  the  field  of  view  of  the  photo¬ 
multiplier  defines  both  the  measurement  volume  and  the  scattering  wave 
vector:  k  *  Hg  -  t^.  The  electric  field  at  the  photocathode  resulting 
from  scattering  by  N  particles  can  be  written 


E(t) 


l  E . [ exp  i(£»r,(t) 
j-1  J  J 


(1) 


where  Ej  is  the  electric  field  amplitude  scattered  into  the  direction 
of  it  by  particle  j,  which  is  at  position  r  .  The  intensity  of  the  scat- 
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tered  light  is  the  ensemble  average  (denoted  by  <  >)  of  the  square  of  the 


electric  field: 


yne  laser-scattering  DBS  experiment. 
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I(t)  -  <E  (t )E(t)> 


The  resulting  photocurrent  is  proportional  to  the  intensity: 


i(t)  -  BE  (t)E(t) 


O 

where  B  is  a  detector  efficiency  (A-cnr-sec/erg). 

Accounting  for  the  finite  relaxation  time  of  the  photocathode  sur¬ 
face,  and  assuming  that  the  temporal  statistics  of  the  light-scattering 

O 

events  are  Gaussian  (ensured  in  dilute  particle  systems),  the  photo- 
current  correlation  function  can  be  related  to  the  electric  field  correla¬ 
tion  function  by  the  relation 

R.(t)  -  <i(t )i(t-hr  )> 

1 


Bel6(x)  +  B2[l2  +  IV0!2] 


where  e  is  the  electron  charge,  6 (t )  is  the  Dirac  delta  function,  and 


R£(t)  -  <E  (t)E(t+r)> 


is  the  electric  field  correlation  function. 

According  to  the  Wiener-Khinchine  theorem,  the  photocurrent  corre¬ 
lation  function  R^(x  )  and  the  photocurrent  power  spectrum  S^(u>)  are  a 
Fourier  transform  pair.  Thus,  from  Eq.  (3), 

1  * 

Si  (a» )  “  —  /  R.  (t  )  exp  (iuit  )  dr 

0) 

-  ±  BEI  +  26  2  I25  (a>  )  +  (B  2/r  )  /  |R£(t  )  | 2  exp  (iurr  )  dr 

—CO 

The  three  terms  on  the  right-hand  side  of  Eqs.  (3)  and  (5)  are  shot  noise, 
the  dc  component,  and  the  ac  term  of  interest  that  depends  on  R^r), 
respectively. 


The  electric  field  correlation  function  for  scattered  light  from  a 
system  of  particles  carried  by  a  fluid  flow  must  still  be  derived.  Using 
Eq.  (1)  in  Eq.  (A)  and  noting  that  the  relative  positions  of  neighboring 
particles  are  uncorrelated  (also  ensured  in  dilute  systems),  one  finds 
that  Rg(x )  depends  on  the  displacement  of  individual  particles  in  time  t : 


VT) 


•  l  <Ea  exp{  ik*  [  r .  (t  )  -  r,(0)] 
j-1  3  J  3 


(6) 


Equation  (6)  is  a  general  expression  for  the  electric  field  correlation 
function  for  scattered  light  from  any  moving  particle.  The  following 
subsections  derive  the  resulting  photocurrent  correlation  function  and 
power  spectrum  for  the  special  case  of  diffusing  particles  in  flowing 
systems.  The  effects  of  high  velocity,  turbulence,  polydispersity ,  and 
short  observation  times  on  particle  size  measurement  are  derived 
independently. 

1.  PURE  DIFFUSION,  POLYDISPERSE  SYSTEMS,  LAMINAR  FLOW 

Displacement  due  to  diffusive  motion  of  the  particles  enables  their 
size  to  be  inferred.  The  probability  P  of  diffusive  displacement  between 
Ar  and  Ar  +  dAr  in  time  x  is  given  by  the  Gaussian  random  walk 


+  -1/2  |Ar  Ct  )  1 2 

d(  A  r  j  (x  ))  -  (AirDjT)  [exp - — ■]  d(Arj(x)) 


(7) 


where  Dj  is  the  diffusion  coefficient  for  particle  j.  The  relation 
between  particle  size  and  the  diffusion  coefficient  is  discussed  below. 
Averaging  the  diffusive  displacements,  one  obtains 

N  9 

x)  -  exp  (-iwgT )  l  Ij  exp[-(|£|  D^  -  (8) 
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where  Ij  is  the  intensity  scattered  by  particle  j  into  the  solid  angle 
collected  by  the  detector.  In  the  absence  of  velocity  gradients  in  the 
scattering  volume,  v^  ■  v^  »  v,  and  the  particle  displacement  due  to 
translation  is  the  same  for  every  particle.  This  Doppler-shif ts  the  opti¬ 
cal  spectrum  by  it  •  v  but  does  not  broaden  the  spectrum  (nor  dampen  the 
correlation  function).  The  resulting  photocurrent  correlation  function 
and  power  spectrum  for  purely  diffusive  motion  are  (neglecting  dc  and  shot 
noise) 


and 


(9) 


(10) 


respectively.  Particle  diffusion  thus  dampens  the  correlation  function  of 
Eq.  (9)  and  broadens  the  power  spectrum  of  Eq.  (10)  according  to  the  val¬ 
ues  of  the  diffusion  coefficients.  A  mean  particle  size  can  be  inferred 
from  the  half-width  of  either  function,  but  as  these  equations  show,  both 
functions  contain  information  about  the  particle  size  distribution  that 
was  responsible  for  scattering  the  laser  light.  Practical  considerations 
are  discussed  below  for  the  inversion  of  scattering  data  to  infer  both  a 
mean  size  and  the  width  of  an  assumed  log-normal  size  distribution. 

The  particle  size  distribution  determines  the  range  of  diffusion 
coefficients  used  in  the  above  sums,  but,  more  Important,  it  determines 
the  weight  given  each  contribution  via  the  scattered  intensities  Ij. 

Since  Ij  is  typically  proportional  to  the  sixth  power  of  the  particle 
radius,  the  small  particles  in  a  wide  distribution  cannot  be  measured. 
Therefore,  the  practical  measurement  of  log-normal  particle  size 
distributions  is  limited  to  those  with  standard  deviations  less  than 
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1.4.^  in  general,  for  polydisperse  systems  without  velocity  gradients, 
the  correlation  function  Is  a  sum  of  decaying  exponentials,  *  and  the 
power  spectrum  is  a  sum  of  Lorentzian  contributions.^’^ 

12  22 

The  method  of  cumulants  »  is  a  convenient  data  reduction  technique 
that  yields  the  coefficient  of  the  first  two  terms  in  the  MacLaurin  expan¬ 
sion  of  the  sum  of  decaying  exponentials.  Comparing  measurements  of  the 
photocurrent  correlation  function  of  light  scattered  from  a  flame  with  Gq. 
(9),  Scrivner  and  coworkers  J  applied  the  method  of  cumulants  to  determine 
both  the  mean  particle  radius  and  the  width  of  an  assumed  log-normal  soot 
particle  size  distribution.  This  method  fits  the  entire  curve  of  the 
measured  correlation  function  to  a  polynomial  expansion  in  t.  In  their 
experiments,  these  authors  used  a  digital  correlator  (commercially  avail- 

Q 

able  with  up  to  5  x  10°  Hz  sampling  rates),  which  efficiently  meas¬ 
ures  R^(t)  as  the  correlation  between  arrival  times  of  photons  at  the 
photocathode.  Note  that  measurements  of  S^Go)  with  an  appropriate  spec¬ 
trum  analyzer  contain  the  same  information  as  measurements  of  R^Cr);  many 
spectrum  analyzers  compute  S^Gii)  by  means  of  fast-Fourier  transformation 
of  the  digitally  correlated  input  signal.  Fitting  the  entire  curve  of  the 
measured  power  spectrum  to  theoretical  curves  based  on  Eq.  (10)  is 
expected  to  yield  accurate  particle  size  distributions.  Using  only  the 
half-width  and  detectability  of  the  measured  power  spectrum  has  been 
shown**  to  be  inadequate  for  measuring  size  distribution,  except  in  a  few 
special  cases.  The  available  hardware  and  software  currently  favor  meas¬ 
uring  the  correlation  function  and  using  the  method  of  cumulants  to  diag- 

A  / 

nose  reasonably  narrow  polydispersity  in  flames. 

2.  RELATION  BETWEEN  DIFFUSION  COEFFICIENT  AND  PARTICLE  RADIUS 

DBS  measures  the  diffusion  coefficient  D  (or  the  distribution  of 
diffusion  coefficients  D ^ )  experimentally.  This  measurement  can  be  theo¬ 
retically  related  to  the  radius  r  (or  distribution  of  radii  rj)  of  the 
diffusing  particles.  The  elementary  diffusion  theory  used  here  assumes 
the  particles  to  be  spherical. 


Types  of  diffusion  are  characterized  by  the  Knudsen  number,  1^,  the 
ratio  of  the  gas  mean  free  path  to  the  particle  radius: 


where  p  is  the  gas  viscosity,  p  is  density,  M  is  molecular  weight,  and  T 
Is  temperature.  For  the  relatively  large  particles  encountered  late  in 


sooting  flames,  Kn  is  small  and  there  is  Stokes -Einstein  diffusion  for 


which 


^Stokes -Einstein 


where  k^  is  Boltzmann's  constant. 

Intermediate  Knudsen  numbers  are  often  encountered  in  flames,  a 
regime  in  which  it  is  appropriate  to  Include  the  experimentally  derived 
Cunningham  slip-correction  factor  in  the  above  relation: 

k  T 

^Cunningham  “  6Wvir  ^  +  +  ^  exP("  iT-)]} 

n 

in  which  A,  Q,  and  b  are  of  the  order  1.2,  0.4,  and  1,  respectively . ^ 
Particles  smaller  than  the  gas  mean  free  path  are  known  to  occur  early  in 
sooting  flames  but  have  seldom  been  measured  by  DBS  because  of  their  small 
scattering  coefficient.  This  regime,  which  has  a  large  Knudsen  number,  is 
known  as  free-molecular,  or  Epstein,  diffusion,  for  which 

k.  TK 

D  -  b  n 

Epstein  4irpr 


3.  PURE  DIFFUSION,  M0N0DISPERSE  AEROSOL,  NO  VELOCITY  GRADIENTS 

Consider  a  monodisperse  aerosol  embedded  in  a  laminar  flow.  The 
single  particle  size  ensures  identical  scattered  intensities  Ij  *  Ik  and 
diffusion  coefficients  Dj  ■  D^.  In  this  case,  the  photocurrent  correla¬ 
tion  function  is  a  simple  exponential  and  the  power  spectrum  is  purely 
Lorentzian:  Eqs.  (9)  and  (10)  become 


R|(t) 
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The  diffusion  coefficient  is  easily  extracted  from  the  half-width  of  the 
measured  function. ^  Note  that  this  measurement  can  be  made  without 


knowing  the  complex  index  of  refraction  of  the  particles,  although  it  must 


be  known  for  size  distribution  measurements  (in  the  calculation  of  Ij  and 


1^  for  the  given  scattering  angle). 


4.  EFFECT  OF  PARTICLE  TRANSIT  TIME  (MONODISPERSE  SYSTEMS) 


If,  because  of  a  high  enough  flow  velocity  v,  the  particle  traverses 
the  measurement  volume  before  it  can  complete  several  cycles  of  its 
Brownian  motion,  the  diffusion  coefficient  cannot  be  measured.  The  par¬ 
ticle  transit  time  is  taken  into  account1^* ^  by  assuming  that  light  is 
scattered  and  collected  from  a  spherically  symmetric  Gaussian  intensity 
distribution  with  characteristic  size  o .  With  the  origin  at  the  center  of 
the  scattering  volume,  the  amplitude  distribution  function  of  the  input 
laser  field  is  given  by 


A(r)-  (2ira2)  exp(-^-l-) 

^2 


so  that  instead  of  Eq.  (1),  the  scattered  electric  field  is 


E(t)  -  Ej[  exp  i{&rj(t)  -  u)Qt}  ]  A(  r^  (t)) 


Averaging  over  both  diffusive  displacement  and  random  particle  position 
gives  an  electric  field  correlation  function  of 
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MO  -  exp(-in)QT  )  i  I  exp{-{  — — +  (  |fc|  D  “itv  )t]) 
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and  modifies  the  photocurrent  correlation  function  of  Eq.  (11)  to  be 

r’(t)  .§i^!l!exp|^Jz|iT2  +  2 1^| 2  Dr]} 
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while  the  power  spectrum  of  Eq.  (12)  becomes  the  Voigt  function^ — the 
convolution  of  Gaussian  and  Lorentzian  functions.  As  discussed  in  Refs. 
15  and  16,  the  particle  transit  time  places  an  upper  limit  on  particle 
sizes  that  can  be  measured  using  a  beam  size  o  in  the  face  of  a  flow  vel¬ 
ocity  v;  particles  have  measurable  diffusion  coefficients  D  if  the  Voigt 
parameter  z  satisfies 

*«(*>  “  ~T  C  2|^|2D)  >  0.2 


In  most  experiments,  tc.  and  a  have  been  made  large  enough  to  require  small 
transit-time  corrections.  Transit-time  contribution  to  the  power  spectrum 
for  polydisperse  systems  is  derived  in  Ref.  16. 

5.  EFFECT  OF  TURBULENCE  (INCLUDING  TRANSIT  TIME,  MONODISPERSE  SYSTEMS) 
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As  noted  by  Chang  and  Penner,  if  the  length  scale  of  turbulence  is 
large  compared  with  the  scattering  volume,  it  has  no  measurable  effect  on 
the  homodyne  power  spectrum  or  correlation  function.  There  is  no  turbu¬ 
lence  effect  because,  at  any  instant,  there  are  no  velocity  gradients  in 
the  scattering  volume,  and  the  results  for  diffusion  in  laminar  flow 
apply. 

If  its  length  scale  is  small,  turbulence  can  produce  a  random  dis¬ 
placement  of  particles  while  they  are  in  the  measuring  volume.  The 
resulting  dampening  of  R^(t)  [broadening  of  S^(u>)]  may  impair  particle 
size  measurements.  Following  Ref.  20,  the  velocity  of  particle  j  is  sep¬ 
arated  into  mean  and  fluctuating  components: 


-  .*  V.V-VJWYiV/W  TV  J*'*  ^ 


♦  +  ■+* 
v.  ■  v  +  Av , 
j  mean  j 


An  isotropic  Gaussian  fluctuating  velocity  distribution  is  assumed,  for 
which  the  probability  of  particle  j  having  a  fluctuating  component 
between  AVj  and  Av^  +  dA v^  is 


P  d(AVjT)  -  < 2n  v2  x2)  3^2  {  exp[ -|AVj  |  2x  2/(  2v2x  2)  ] }  d(AVjX  ) 
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where  Vj  is  the  characteristic  width  of  the  velocity  fluctuation  distribu¬ 
tion.  In  calculating  the  electric  field  correlation  function  the  ensemble 
average  over  turbulent  displacement  Av  is  included  with  the  average  over 
diffusive  displacement  r,  and  the  average  over  particle  position  r  within 
the  volume  defined  by  A(r).  The  result  is^® 


2  2  2 
v„x 
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RjXx  )  *  exp[  i(£*  v-u0)x]  l  I^[  exp-(  |^|2DjX  +  - ^ ^ ^ 


The  photocurrent  correlation  function  becomes 

RJ(t)  -£-~-  exp{-[2|^j2Dr  +  (|£|2vT2  + 

Ao 

and  the  photocurrent  power  spectrum  is  "the  convolution  of  the  following 

three  components:  a  Lorentzian  profile  due  to  the  Brownian  motion,  a 

Gaussian  profile  due  to  the  turbulent  velocity  fluctuations,  and  another 

Gaussian  profile  due  to  the  nonuniformity  of  illumination  in  the  coherence 
20 

volume."  This  spectrum  is  again  a  Voigt  function  whose  Lorentzian  part 
2 

is  still  2K  D,  but  its  Gaussian  component  is  now  proportional  to 


[4|£|2  v2  +  (v2/o2)]~1/2 


The  Gaussian  component  creates  an  upper  limit  of  particle  sizes  that  can 
be  measured  in  the  face  of  turbulence,  similarly  to  the  upper  limit 
Imposed  by  transit  time  discussed  in  the  preceding  subsection.  Mono- 
disperse  particles  can  be  sized  if  the  experimental  conditions  are  such 
that 


Re(z) 
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which  results  in  a  measurable  Lorentzian  contribution.  Chang  and  Penner 
note  that,  for  certain  combinations  of  particle  size,  measurement  volume, 
and  turbulence  level,  three  separate  homodyne  measurements  at  different 
scattering  angles  can  yield  v,  v^,  and  D  simultaneously.  Heterodyne  LDV 
measurements  of  v  and  vT  can  always  be  used  to  accurately  assess  the 
Gaussian  contributions  and  to  determine  the  feasibility  of  particle  size 
measurement. 


6.  EFFECT  OF  SHORT  OBSERVATION  TIMES 

To  use  DBS  for  time-dependent  particle  size  measurement,  or  for  meas¬ 
urement  using  pulsed  laser  sources,  the  effect  of  short  observation  times 

must  be  investigated.  This  effect  is  discussed  in  detail  by  Penner, 
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Bernard,  and  Chang.  The  simplest  analysis  assumes  that  the  observation 
time  is  much  shorter  than  the  particle  transit  time,  and  that  turbulence 
can  be  neglected.  Here,  the  scattered  electric  field  is  multiplied  by  a 
square  function  II  (t)  representing  the  observation  time  T: 


II(t;T)  - 


1  for  0  <  t  <  T 
1/2  for  t  -  0  and  T 
0  for  t  >  T 
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The  autocorrelation  of  the  square  function  is  the  triangle  function 


1  -  t/T  for  0  <  t  <  T 
A(t;T)  *  1  for  t  *  0 

0  for  t  >  T 

so  that  the  scattered  field  correlation  function  for  a  single  particle 
size  becomes 

Rj.(t)  -  NI  A(t;T)  exp(-imQT)  exp(-|ic|2Dr) 

The  photocurrent  correlation  function  is  damped  by  the  triangle  function 


squared: 


Rt(T) 
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and  its  power  spectrum  becomes  the  convolution  of  the  Lorentzian  diffusion 

)  6 

contribution  and  the  spectrum  of  the  square  of  the  triangle  function. 

The  latter  contribution  is 
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Experimenters  must  determine  which  of  the  corrections  discussed  in 
Subsections  4  through  6  above  apply  to  their  DBS  measurement  (discussed  in 
Subsections  1  and  3),  and  minimize  those  corrections  over  which  they  have 
control. 


B.  SELECTED  EXPERIMENTAL  RESULTS 

DBS  has  been  used  to  measure  soot  particle  sizes  in  many  types  of 
flames.  Studies  on  flat  flames  of  ethylene,  acetylene,  methane,  and  pro¬ 
pane,  and  on  wick  flames  of  kerosene  and  paraffin,  appear  in  published 
literature.^  Two  recent  studies are  worthy  of  review  here  because 
they  demonstrate  the  maturity  of  this  technique  and  its  utility  in  eluci¬ 
dating  the  mechanism  of  soot  growth. 
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Weill,  Lhulssier,  and  Gouesbet^  have  measured  mean  particle  radii 
above  a  methane-oxygen  flat  flame  using  DBS,  and  have  also  Inferred  soot 
number  densities  at  each  measurement  location,  based  on  the  scattered 
signal  strength  and  the  measured  particle  size.  These  data  were  taken  as 
a  function  of  height  above  the  burner  for  two  fuel/oxygen  equivalence 
ratios,  with  three  values  of  the  unburned  gas  velocity  at  each  equivalence 
ratio.  An  interesting  result  is  obtained  when  the  particle  size  and  num¬ 
ber  density  are  plotted  as  a  function  of  residence  time  in  the  flame, 
Instead  of  height  above  the  burner:  Regardless  of  the  unburned  gas  veloc¬ 
ity,  flames  of  the  same  equivalence  ratio  produce  identical  soots  at  iden¬ 
tical  times  in  the  flame. This  result  is  expected  for  a  soot  growth 
mechanism  (such  as  surface  growth  or  coagulation)  in  which  the  particle 
size  depends  on  the  number  of  collisions  it  has  encountered. 
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Particle  size  and  number  density  can  be  used  to  calculate  the  soot 
volume  fraction  at  the  measurement  location.  Weill  et  al.  have  noted  that 
if  surface  growth  were  responsible  for  the  increase  in  particle  size  with 
height  above  the  burner,  the  number  density  would  remain  constant  as  the 
soot  volume  fraction  increased,  and  coagulation  of  smaller  particles  would 
result  in  decreasing  number  density  and  constant  volume  fraction.  Their 
plot  of  number  density  versus  soot  volume  fraction  is  depicted  in  Fig. 

2.  For  the  two  flames  with  higher  flow  velocity,  a  nearly  constant  volume 
fraction  accompanies  the  decrease  in  number  density  with  height  above  the 
burner,  suggesting  the  dominance  of  coagulation  for  soot  growth  in  these 
flames.  Only  very  late  in  the  faster  flames,  or  throughout  the  slow 
flame,  when  the  number  density  is  well  below  10+^  m”^,  does  surface 
growth  become  evident. 

2S 

Scrivner,  Taylor,  Sorensen,  and  Merklin*"'  have  measured  both  the  mean 
particle  radius  and  the  width  of  the  soot  particle-size  distribution  as 
functions  of  height  above  a  premixed  methane /oxygen  flat  flame.  They  fit 
the  measured  correlation  function  to  a  second  order  expansion  in  t: 


*w 
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R^(t  )*  Cq  exp(-2K^T  +  K2t2) 

The  mean  pardcle  radius  and  the  standard  deviation  of  the  assumed  log¬ 
normal  particle  size  distribution  were  obtained  from  Kj  and  K2,  the  first 
two  cumulants  in  this  expansion.  Scrlvner  et  al.'s  measurements  show  a 
decrease  in  the  standard  deviation  of  the  distribution  with  height  above 
the  burner,  which  they  interpreted  as  verifying  the  dominance  of  coagula¬ 
tion  in  soot  particle  growth.  Finally,  they  quote  the  theory  of  diffu- 
sion-controlled  coagulation  (developed  by  Pich,  Friedlander,  and  Lai4  ), 
which  predicts  the  narrowing  of  an  initially  broad  distribution  toward  the 
limit  of  the  approximately  log-normal  self-preserving  size  distribution. 

This  corroborates  the  light-scattering  and  DBS  measurements  of  Chang  and 
20 

Penner  that  also  gave  evidence  that  the  self-preserving  distribution  is 
attained  in  sooting  flames. 
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III.  IMAGING  TECHNIQUES 


Particles  larger  than  a  few  microns  are  difficult  to  size  by  DBS,  but 
can  be  accurately  sized  by  microscopy  or  holography.  Motion  pictures  can 
be  made  of  droplet  and  particle  combustion,  and  large  soots  from  combus¬ 
tion  exhaust  can  also  be  photographed.  The  complicated  optical  systems 
required  for  microscopy,  and  their  small  depth  of  field,  can  be  overcome 
by  holography.  Holographic  measurement  of  particle  size,  shape,  and 
velocity  is  discussed  in  Section  III. A.  In  Section  III.B  are  described 
techniques  by  which  transparent  droplets  can  be  made  to  emit  laser  light 
or  Raman-scattered  light  to  enable  them  to  be  photographed  easily. 

A.  HOLOGRAPHIC  PARTICLE  DIAGNOSTICS 

Like  DBS,  far-field  holography  was  developed  soon  after  the  laser  was 
invented.  Its  first  description  and  considerable  development  can  be 
attributed  to  B.  J.  Thompson.  Far-field  holography  has  become  an  accepted 
technique  in  particle-size  analysis  and  is  complementary  to  DBS:  it  can 
measure  particles  of  1-pm  diameter  and  larger.  It  can  also  yield  the 
particle  velocity  and  the  cross-sectional  shape  of  the  particle. 

Although  like  DBS  in  its  conceptual  simplicity,  such  holography  involves 
data  reduction  that  requires  careful  analysis  of  the  optical  and  flow 
characteristics  of  the  experiment. 

A  schematic  diagram  of  the  recording  of  a  far-field  hologram  is  shown 
in  Fig.  3a.  A  photographic  emulsion  is  illuminated  by  two  samples  of  a 
short  pulse  of  laser  light.  The  first  sample  is  a  plane  wave  (a  beam 
whose  surfaces  of  uniform  phase  are  flat)  taken  from  the  expanded  laser 
beam,  and  the  second  sample  is  light  diffracted  from  the  same  beam  by 
particles  within  the  measurement  volume.  The  interference  of  these  two 
beams  in  the  plane  of  the  emulsion  is  recorded  as  a  hologram.  The  holo¬ 
gram  is  reconstructed,  as  depicted  in  Fig.  3b,  by  being  illuminated  with  a 
plane  wave  from  a  continuous-wave  laser.  The  images  of  the  scattering 


particles  can  be  studied  from  the  reconstructed  hologram  or,  alter¬ 
natively,  the  particles  can  be  sized  by  direct  examination  of  the  inter¬ 
ference  patterns  on  the  hologram.  The  practical  aspects  of  this  method 
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are  reviewed  by  Cartwright,  Dunn,  and  Thompson,  including  the  required 
film  resolution  and  format,  laser  power,  coherence,  and  pulse  duration, 
and  the  extent  to  which  particle  shape  can  be  inferred. 

Exposing  the  hologram  with  two  short  pulses  enables  the  particle 
displacement  between  pulses  to  be  measured  and  hence  the  particle's  aver¬ 
age  velocity  during  that  time.  A  novel  optical  data  processing  technique 

for  such  double-exposure  holograms  is  described  by  Malyak  and 
10 

Thompson.  That  technique  directly  yields  the  velocity  distribution 
within  the  measurement  volume. 

B.  RAMAN  SCATTERING  AND  LASING  DROPLETS 

After  a  light  field  has  been  applied  to  them,  transparent  droplets  or 
particles  (or  fibers)  can  exhibit  internal  resonances,  depending  on  the 
particle  size  and  refractive  index,  and  the  light  wavelength.  When  such 
resonances  occur  when  internally  reflected  light  travels  an  integral  num¬ 
ber  of  wavelengths  during  one  round  trip  inside  the  droplet,  they  are 
known  as  morphology-dependent  resonances.  Counterpropagating  resonant 
waves  can  result  in  a  standing  wave  inside  the  droplet  with  a  high  inten¬ 
sity  near  its  surface.  The  droplet  surface  appears  to  be  highlighted  by 
light  scattered  from  this  wave,  and  photographs  of  the  droplet  exhibit 
clear  outlines  of  its  shape.  Studied  extensively  by  R.  K.  Chang  and 
coworkers,  the  light  emitted  by  the  droplet  can  be  elastic  scattering, 
spontaneous  Raman  scattering,  fluorescence,  stimulated  Raman  scattering 
(up  to  the  fourteenth  order),  or  laser  emission.  For  laser  emission  to 
occur,  the  wave  must  experience  a  net  gain  during  one  round  trip  of  the 
droplet.  Such  a  gain  was  achieved  by  Chang  and  coworkers  by  illuminating 
rhodamine-doped  ethanol  droplets  with  pulses  of  532  nm  wavelength  and 
observing  the  dye-laser  output  at  595  to  600  nm.  Notable  from  their 


photographs  is  that  the  droplets  need  not  be  spherical  or  of  a  particular 
size  to  exhibit  laser  emission.  Indeed,  their  photographs  of  highly  dis¬ 
torted  droplets  and  liquid  streams  demonstrate  surface  lasing  to  be  a 
powerful  technique  for  highlighting  droplets  to  study  their  dynamics. 


IV.  CONCLUSION 


Laser  diagnostics  of  particle  size  have  been  reviewed,  emphasizing 
their  use  in  combustion  systems.  For  the  small  particles  characteristic 
of  soot  in  flames,  the  technique  of  diffusion-broadening  spectroscopy  was 
discussed  in  detail.  This  in  situ  technique  has  also  proven  to  be  applic¬ 
able  to  measuring  distributions  of  particle  sizes  for  which  the  standard 
deviation  is  less  than  1.45.  It  has  yielded  the  important  conclusion  that 
soot  particles  coagulate  into  a  self-preserving  size  distribution  within 
the  flame.  Two  distinctly  different  laser-based  techniques  for  Imaging 
larger  particles  were  described.  These  techniques  have  not  yet  been 
applied  to  combustion  systems,  but  such  an  application  Is  certainly  feas¬ 
ible.  Further  discussion  of  the  topics  covered  here  is  available  in  the 
cited  references  and  in  their  bibliographies. 
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